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NATIOML ADVISORY COMMIMEE FOR AERONAUTICS 


RESEARCH MSMDRANDUM 


SOME LOW-SiEED STUDIES OF THE EFFECTS OF WING LOCATION 
ON WING-DEFORMATION— BODY-FREEDOM FLUTTER 
By E. Widmayer, Jr. 


SUMMAET 


An investigation of flutter of wings moxinted on the after portion of 
a body free to pitch was made in the Langley 4. 5*foot flutter research 
tunnel. The purpose of the investigation was to explore the effect of 
wing location on a low-frequency flutter which had occurred in some rocket 
tests. Uniform wings of 0^ and 45^ sweep were tested at various positions 
rearward of the pitch axis. All experiments were performed in air at 
atmospheric pressure and at low speeds. 

A low-frequency flutter involving primarily wing deformation and "body 
pitching was obtained on the unswept wings. The flutter frequencies were 
in the order of the ** short period** stability frequencies, being consid- 
erably below the lowest natural deformation frequencies of the wing. The 
reduced frequencies at flutter were in range of 0.01 to 0.03, representing 
a wave length of from 210 to 500 semichords. No low-frequency flutter was 
obtained with the 45^ swept wings for the few cases considered. 

A comparison of the experiment has been made with theory, with various 
assumptions regarding the air forces being used. The experimental trend 
of the flutter-speed coefficient to Increase as the unswept wings were 
moved rearward was also obtained from the analysis. It was found that 
taking account of finite span effects gave better agreement between exper- 
iment and theory than did the use of two-dimensional air forces. The 
experimental and calculated frequencies of flutter -vdxen referred to the 
natural deformation frequencies of the wing were of the same order of 
magnitude . 


INTRODUCTION 


During some previous flutter tests made with rocket models having the 
wings rearward of the configuration center of gravity, a low-frequency 
flutter involving wing deformation and . rigid-body freedoms was encountered 
on unswept wings (ref. 1). An analysis of this type of flutter indicates 
that the mode of flutter is strongly influenced by the wing location with 
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respect to the configuration center of gravity. The problem of the effect 
of rigid-body freedoms (fuselage mobility) on flutter has been of Interest 
for some time (refs. 2 to 4 ). The early Investigations for the most part 
dealt with conventional locations of unswept wings, the system center of 
gravity being close to the wing quarter chord. Interest in the subject 
has increased with the use of swept wings, tailless configurations, and 
wing locations remote from the system center of gravity (refs. 5 to 8). 
References 6, 7^ and 8 treat experimentally flutter of swept, cranked, 
and delta wings having root freedoms in pitch and translation. 

In some recent wind-tunnel experiments on the flutter of wing and 
horizontal-tall combinations, a low-frequency flutter involving wing 
deformation and body pitching was observed when the tail was absent. 

Since flutter of this kind may have a bearing on the stability of missiles 
and tailless designs and since there seems to be a lack of data on this 
type of flutter, these experiments on a tailless model are being reported. 
These experiments, made at low speeds, were performed for some uniform 
wings having 0 ° and 45 ° sweep in which the wing was located rearward of 
the pitching axis. A comparison is made of all results of the experi- 
ments with some analytical results in which aerodynamic forces based on 
various approximations were used. The flutter involved low reduced fre- 
quencies, and it might be expected, that the effects of aspect ratio would 
be important. The calculations were based on the following air forces; 

(a) the unsteady two-dimensional theory of reference (b) the \msteady 
two-dimensional theory of reference 9 modified to include finite-span 
effects as derived in reference 10, (c) the air forces of reference 9 for 
the case of C(k) =1, (d) the air forces of (c) modified by the correc- 
tion a/(a + 2) . 


SYMBOIS 


A aerodynamic aspect ratio, for rectangular plan- form wings, Z/c 

a nondlmensional distance (in wing semichords) from wing midchord 

line to location of the elastic axis, positive when elastic 
axis is behind wing midchord line 

I wing span 

C(k) Theodorsen’s function F(k) + iG(k) (ref. 9 ) 

c wing chord 

dfj expressions appearing ;Ln determinant (defined in appendix) 
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' first iDending natural frequency^ cps 

f^ first torsion naturaJL frequency, cps 

damping coefficients as defined in reference 11 


h 




k 

m 

0 

q 


s 


V 

2V/c%c 

y 

a ‘ 


the bending deflection of wing at span station y given in 
terms of reference deflection (at tip) and norma^lized mode, 

mass moment of inertia of fuselage and wings about the pitch 
axis 

mass moment of inertia of wing per unit length about wing 
elastic axis 

reduced frequency parameter, ox:/2V 

mass of wing per unit length 

angle of pitch of fuselage center line, positive when nose up 

dynamic pressure 

nondimens ional radius of gyration of wing section about wing 

p ^ w 

section elastic axis, r = 4 — ^ 

nondimens ional distance (in wing semichords) from wing midchord 
line to fuselage pitch axis measuxed parallel to fuselage center 
line, positive when pitch axis is behind wing midchord line 

flutter speed 

nondimens ional flutter speed coefficient 

nondimens ional position of wing-section center of gravity in 
wing semichords relative to the wing section elastic axis, 
positive when center of gravity is behind elastic axis 

coordinate of span station measured from fuselage center line 

torsional deflection of wing at span station y in terms of 
reference deflection a* (at tip) and normalized mode. 
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ki. ratio of mass of wing to the mass of a cylinder of air of a 

diameter equal to the chord of wing, both taken for equal 

length along the span, n = — 

jtpc^ 


P 

°CL» 

a 

O) 

“h 


mass density per unit volume of air 

correction factors to C(k) which vary with span station as 
given In reference 10 

parameter defined by fl = (coj^/co) (l + ig) 

circular frequency of flutter 

circular frequency of natviral bending deformation 
circular frequency of natural torsional deformation 


APPARATUS AND METHOD 


Tunnel .- The experiments were conducted in air at atmospheric pres- 
sure in the Langley h, 5-foot flutter research tunnel which is a closed- 
throat single-return type liaving a cylindrical test section. The highest 
velocity reported in these tests is 338 feet per second, corresponding 
to a Mach number of 0.32. The Reynolds number was greater than 1 x lOo. 

Model . - The model used in these experiments was intentionally sim- 
plified and is shown in figure 1. The model consisted of a fuselage with 
vertical tail and the aerodynamic surfaces to be examined. The fuselage, 
of 1- inch-square aluminum stock, was 31 .75 Inches in length. The fuse- 
leige was mounted on a ball-bearing pivot and had no spring restraint in the 
pitching degree of freedom. Provision was made for this bearing to be 
located at various positions on the forward 13 Inches of the fuselage. 

Lead counter weights were enqiloyed to maintain, as near as possible, 
static balance of the entire model about the bearing support. No pro- 
vision was made for maintaining a constant mass moment of Inertia about 
the pitch axis. Consequently, this inertia varied with the position of 
the wings and with the counter weighting used for each pitch axis. The 
wings could be attached to the fuselage at six different stations. 

The model was mounted in the t unne l by use of taut wires connected 
to the stationary bearing support^ thus, although the bearing housing 
was held firm, the fuselage could pitch about the bearing axis and had 
some freedom in yaw. A rigid (as compared with the wings) vertical tail 
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of 2.8 X 9-5 X 0.125 aJLuinlnum plate effectively restrained the freedom 
of yaw when rnider the influence of the air stream. 

Wings . - The different sets of wings were mounted in "base blocks to 
facilitate the variation in their location relative to the pitching axis. 
Experiments were conducted on two sets of wood wings, designated wings I 
and wings II. These wings had an airfoil section of NACA 16 -OO 6 and were 
made of solid maple. It was necessary to cut slots perpendicular to the 
plane of the wing in the chordwise direction to reduce the struct-ural 
stiffnesses so that the desired flutter might be obtained at the desired 
speeds. The structural properties are presented in table 1(a). 

The 45 ^ swept wings used in the experiment were ^ ~ inch aluminum- 

alloy sheet. These wings had the same tlp-to-tip distance as wings I 
and II and had a sheared tip plan. The structural properties are pre- 
sented in table l(b). 

Instiumentation . - The instrumentation consisted of resistance-type 
strain gages mounted at the root of each wing so as to indicate the bending 
or torsion deformation or both- The strain-gage output was recorded on 
an oscillograph. The action of both the fuselage and the wings could be 
qualitatively determined from this record, since from observation it was 
noted that the pitching of the fuselage was accompanied by the bending of 
the wings. Although this instrumentation yields frequencies of oscilla- 
tion, the phase between bending and torsion, and relative an^lltudes of 
the wing defoimations, it does not give the phase between the pitching 
of the fuselage and the wing deformation or the amplitude of the pitching 
oscillation. 

Test procedure .- In order to obtain the flutter of the model for a 
given wing location, the airspeed was increased slowly until flutter was 
observed. The strain-gage outputs and the tunnel conditions were then 
noted . 


ANALYTICAL CONSIHE^IATIONS 


The analytical procedures are indicated in this section in order to 
provide a background before presenting the results of the experiments. 

In previous investigations, involving cases where t^e values of reduced 
frequency are fairly high, it has been found that the flutter calculations 
based on two-dimensional oscillating air forces usTially give good agree- 
ment with experiment. In these tests, however, the values of the reduced 
frequency parameter were low, bordering on those encountered in certain 
stability phenomena. It was felt, therefore, that the two-dimensional 
theory would have to be modified to account for aspect-ratio effects and 
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possllDly other factors. Accordingly, compixtatlons of the flutter vere 
made using several different concepts of the air forces and moments 
involved. These concepts are designated as follows: 

A air forces of two d i m ensional incompressible flow as in 

reference 9 

B air forces of concept A modified by the aspect ratio correction 

of Relssner as in reference 10 

C air forces of concept A with the circulation function C(k) = 1 

D air forces of concept C reduced by the aspect ratio correction 

factor 

A + 2 

In A to D only the aerodynamic terms in the equations of equilibrium are 
altered, the structural terms remaining unchanged. 

The model was treated as having three degrees of freedom, pitching 
about the bearing pivot, and two wing-deformation degrees of freedom which 
are approximated by the first uncoupled bending and first uncoupled 
torsion modal shapes of an ideal uniform cantilever beam. In order to 
determine if the structural representation was adeq,uate, a differential- 
equation type of flutter analysis as presented in reference 12 was per- 
formed on one case. The results of this analysis are shown in figure 2 

and indicate that the representation was adequate. 

The wings were treated as though they extended through the fuselage, 
and no attempt was made to approximate the air forces arising from the 
fuselage. The aerodynamic forces at each span station were assumed to be 
proportional to the displacement and motion of the wing section at that 
station (strip analysis). 

Since the pitching degree of freedom has zero natural frequency, the 
equations of motion for the model undergoing harmonic motion are of the 
following form: 



+ di20o + di36o 


= 0 





^22 - 


+ d230o = ° 




2h 

dsi + d3gaQ + d 330 Q = 0 


( 1 ) 
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The equations and the definitions of the quantities dj^j are given in the 
appendix. The “border condition of stability or flutter is given for the 
determinant of the coefficients equal to zero. 


RESULTS AND DISCUSSION 


The results of the experiments and of the computations are presented 
in table II and Illustrated in figures 2 and 3. Computations A, C^ and D 
(see appendix) were performed for each experimental point for unswept 
wings involving the low-frequency wing-deformation — body-pitching flutter 
while computation B was perfomed for two cases. In the experiments using 
45^ swept wings no low-frequency flutter was obtained and no calculations 
were performed. The conditions for these experiments are given in 
table III. 

Unswept wings .- In table II it may be seen that the experimental 
value of 2V/ccD decreases as the distance from the root midchord to 
pitch axis is decreased. A similar trend is also noted for all the 
analytical results. It is seen that the order of magnitude of the 
anal^'tical coefficients for B and D are in agreement with the experi- 
mental values; whereas^ considerable differences are noted between the 
results of A and C and the experimental 2 V/ccdjj coefficients. The 
significance of correcting the air forces for some effects of aspect 
ratio may be evidenced by the improvement in the results with the experi- 
ment when some correction is made. 

It is of interest to note that calculation A indicated that flutter 
of this type was unobtainable for s =s -2.4 while calculations C and D 
predicted flutter. As indicated in figures 2 and 3 the nature of the 
flutter as discerned from analysis is changing to a different type of 
Instability. An attempt was made to obtain wing deformation — body- 
pitching flutter experimentally when s lay in the region described as 
single degree of freedom instability^ but no instability was noted (see 
table II). 

In order to determine the bending-tors ion-flutter characteristics of 
wings some tests were made with the fuselage restrained in pitch. An 
instability of destructive violence was encountered which resulted in the 
loss of the wings. The loss of the wings came without warning and^ iinfor- 
tunately^ only the air velocity at destruction was obtained. From cal- 
culations for wing bending- tors ion flutter using the forces of reference 9 
a theoretical flutter speed was obtained that was 17.4 percent lower than 
the experimental velocity. The torsional divergence speedy as calculated 
by the method of reference 11, is 9 percent lower than the experimental 
velocity. These data are given in table II and figure 2. No tests of 
this kind were made for wings II for they had been damaged in previous 
testing. 
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The correlation of the calculated frequencies to the esgoerimental 
frequencies may be seen by referring to table II. Calculated frequencies 
of Cy and D appear to be consistently hi^er than those of the experi- 
ment, while the calculated frequencies of A appear to be consistently 
lower than the experimental values. However, the trend of the experi- 
mental frequency to increase as the wings are moved rearward is predicted 
by the results of the analysis. Furthermore, if both experimental and 
computed frequencies are referred to the natiiral defonnation frequencies 
of the wing (fj^ = 26.9 cps, fa; = 106 cps), it is seen that the cal- 
culated frequencies are of the same order as the experimental frequencies. 

Swept wings .- In the experiments involving the 45^ swept wings, no 
low-frequency flutter was obtained for the range of wing locations 
examined. The wlng-root-mldchord location was varied from s = 0 to 
s = 5-88, where the reference semlchord is taken normal to the wing 
leading edge. The data for these experiments are presented in table III. 
In references 6, T, and 8 wing flutter Involving body freedoms in both 
pitch and tianslatlon has been reported for the swept wing. In those 
tests, conditions of small body inertias and the mass center location 
within the root chord of the wing were used. In the present sinple 
experiments, these conditions were not obtained and could possibly 
account for the absence of flutter. 


CONCLUDING EEWAEKS 


A low-frequency flutter, involving interaction between rigid-body 
pitching and wing deformation was obtained experimentally on an unswept 
wing. The flutter-speed coefficient increased as the wing was moved 
rearward along the fuselage away from the pitch axis. This trend was 
also indicated by theoretical calculations based on various assumptions 
regarding the air forces. No flutter was obtained for a similar wing- 
body combination with 45^ swept wings for the few cases considered. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Langley Field, Va. 
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APPENDIX 

DETERMINATION OF THE FLUTTER DETERMINANT 


The equations of motion for a system having a body freedom in pitch 
and wing deformation freedoms in twisting and bending are obtained from 
Lagrange's equations. The kinetic energy of the system is as follows: 


1 

T = I IpS 


r ^ Iw(a^ + 2ae)dy 
^0 



2/2 



+ 2(-s + a + h0 + 2xjx ^ ho, + 2(-s 


+ 



dy 


The potential energy is given 

t>y 



of 

'i/2 

u = 0^2 r 
Jo 

mh^dy + J 

IvCt^dy 

) 

The generalized work is given 

by 


-3 o 

r fo-h 

“1 

dh = -2itp ^ j 

% 

r ^ Aeh + 
0 ^ 

eA^e + oAcJcih dy 



“j 

o 

0^ da = -2itp — ^ cs:r 
lo 


0Aae + dy 

c^ o 

r Fib 

“1 

Qfl d9 = -2np — 

^ 1^0 

^0 1!"-^ 

SAjug + d0 dy 
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From Lagrange ' s equation the following equations are obtained: 

> 2 h^ 


(‘^11 - ‘^12“'o + ‘^13® o = 


di 2 

^ + (d 22 - 

0 

tl 

0 

CD 

on 

+ 

X 

di3 

2 h 0 

— + d32ao + d33 

0 

II 

0 

CD 

In these equations ! 

^ ^ 

where g has the 

damping coefficient. 

The quantities 

dj^j are 

dll = 

H - Ach 


^12 = 

(tix^ - Aca)li 


•3-13 = 

[^(-s + a + xix) - 

Ac^ I 2 


d21 = - Aah)(^) I3 

<^2 = 

<^23 = + (-S + a)x, 


=3 - 


14 


"31 


^32 


= [^(-s + a + Xa) - I 5 


+ (-S + a)x, 


= {^E«‘ 

33 - ^ -a.) 


“ -^da f ^ 


2 "6 


where lirg^ = moment of fuselage about pitch axis 


(span)jtp(c/2) ^ 
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The In are, for modal shapes considered. 


So 

f 

lo 


3558 


^h^e <3y 


^3 = 


r 

'Jq 

r 

n 


= 1.5662 


%^<3y 


faf h <3y 

^ 0.67780 


/ 

’-'0 


f a dy 


It. = 


Ve <3y 

^0 


r 


n 


= 1.T2T3 


fcT<3y 




l6 = 


/' 

^0 

'70 


= 0.39153 


fe^dy 




r 

^0 


= 0.63662 


f0^<ay 


The Aj^j are the aerodynamic expressions which are modified in 
various ways as depicted in the section entitled "Analytical 
Considerations" and are given as follows: 


Analysis A .- 
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^c9 ~ 

^^0 ~ ~ ^act “ (s - a)Agj2 

^dh = Ag^(l^^s) 

Ada = A^(k,a,s) = A^^^ - (s - 
Ad0 = Aaa(k,s) 

Analysis B .- 
Ach = Ach + 

Aca “ A^a ■*" "*" ®) 

^ = Agj, - a^{y)^(^ + a) 

Aaa = Aao, + 0a(y) + a^ - l^i - e^ 

^ = Aee + «^e(y)jj ^ - ®)] 

A^ = A^ + ffjj(y)^ + s) 

^ = Aae - '’e(y)j^(| -^ ^) + f-g - =^)| 

^ = Aaa - “e(y)j^§ + •) + - "-2^ 

I 

(s - a)0Q(y)^^ + a) 


- (s - a)Agj^ + 
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Ado, - Aaa - '’o(y)|^^ * a) + ^(1 - s2)j - (s - a)Aco - 
(a - a)da(y)[^ + ¥-(| - ^ 

Analysis C »- Evaltiated as for analysis A for F = 1, G = 0. 
Analysis D.- of analysis C. 
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TABLE I.- WING PARAMETEIS 
(a) Unswept Wing 




Wing I 

— 
Wing II 



Ri^t 

Left 

Right 

Left 


Section 

NACA 16-006 

RACA 16-006 

NACA 16-006 

NACA 16-006 


Material 

itaple 

Maple 

Maple 

Maple 

Semi span length (Center line 
of fuselage to tip) 

1.298 ft 

1.290 ft 

1.298 ft 

1.298 ft 


Density 

0.0216 Ih/in.^ 

0.0216 l-b/ln .3 

0.0216 lb/ln .3 

0.02i 6 lij/iii .3 


c 

0.333 ft 

0.333 ft 

0.333 ft 

0.333 ft 


Elastic axis 

^ 3-35 percent c 

38.0 percent c 

45 . 1 *’ percent c 

40,6 percent c 

-Section center of gravity 

l^■ 6-8 percent c 

1 * 6.8 percent c 

46.8 percent c 

46.8 percent c 


a 

-0.133 

- 0 . 2 l* 

-0.092 

-0.188 


% 

0.069 

o.it 6 

0.028 

0.124 



0 . 23 ^^ 

0.261 

0.230 

0.245 


g 

0.03 

0.03 

0,03 

0.03 


fh 

29.9 cps 

29.9 cps 

27.0 cps 

27.0 cps 


fa 

106 cps 

106 cps 

105 cps 

105 cps 


” j A 

7.25 

7.25 






1l3I3Sl m ■tfOYH 
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TABLE I.- WING PARAMEEERS - Concluded 
(Td) 45 ° Swept Wings 


Section 0 . 06 i(— Inch- thick flat plate 


Material alumln\m 

Semispan length (center line of fuselage to tip) ft 1.298 

Chord perpendicular to leading edge, ft 0.333 

Density, Ih/ln.^ 0.101 

Elastic ELxls, percent chord ^ 

Center of gravity, percent chord 50 

f^, cps 6.1 


fa, cps 


33.3 



TABIJB !!•- FIOTTER DA!I!A 


•» 


ik 


CO 


X 


(a) WlngB 1, fuBelage free 


V 

ft/sec 

q 

Ib/ft® 

I 

s 

2V/ood 

Experimental 

2V/aoh 

Sxperlaantal 

Analytical 2V/oo3i^ 

^CpB 

Analytical f{,ps 

A 

B 

c 

D 

Experimental 

A 

B 

c 

D 

2>(8.1 

73.7 

0.2725 

-2.412 

79.60 

8.86 

(a) 

(D) 

6.85 

^ 7.72 

2.98 

(a) 

(b) 

3.09 

3.08 

259.3 

30.6 

.1800 

-4.68 

42.20 

9.21 

5.17 

(b) 

7.85 

8.84 

5.88 

3.95 

(b) 

6.91 

6.88 

282.1 

95.1 

.2180 

- 6.43 

40,65 

10.09 

6.55 

8.95 

8.4o 

9.45 

6.64 

5.67 

9.5 

8.38 

8.41 

282.3 

95.1 

.2950 

-7-50 

39.00 

10.10 . 

7.30 

(b) 

8.67 

9.90 

6.93 

6.77 

(b) 

6.20 

9.17 

292.6 

101e9 

- .2650 

- 8.26 

34.25 

10.45 

7.34 

(b) 

8.38 

9.90 

8.16 

6.89 

(bi 

9.71 

9.60 

2^.6 

99.8 

.2790 

- 8.26 

33.45 

10.34 

7.30 

(b) 

8.81 

30.30 

8.28 

6.82 

{b 

9.39 

10.22 

297.1 

104.7 

.3067 

—3.26 

39.00 

10.61 

7.26 

w 

8.70 

9.74 

7.28 

6.22 

(b) 

8.66 

9.06 

295.6 

103.8 

.3560 

- 8.26 

39.70 

10,56 

7*20 

(b) 

B .60 

9.87 

7 b 12 

5.92 

(b) 

8.31 

8.23 

297.7 

105.2 

.3425 

- 9.32 

38.80 

10.63 

7.65 

9.65 

9.03 

10.15 

7.33 

6.97 

10.2 

9.55 

9.59 






Tusalage Fixed 









338.1 

137.1 

(a) 

(d) 

(c) 

12.08 

10.00 

(D) 

(b) 

(b) 

(c) 

56.4 

(b) 

(b) 

(b) 

(b) Vingfl II , fuselage free 

283.1 

93.45 

0.4513 

- 11.25 

36.40 

10.19 

8.15 

(b] 


9.35 

10.48 

7.49 

7.54 

(b) 

9.88 

9.83 

271.7 

85.8 

.3522 

-9.42 

35.05 

9.71 

7.72 

Cb 


9.1 

10.15 

7-.41 

7.16 

(b) 

9.68 

9.43 

259.0 

77.5 

.2927 

- 7.65 

38.65 

9.26 

7.05 

w 


8.68 

9.75 

6.41 

6.02 

(bj 

8.4l 

8.57 

226.7 

57.4 

.2041 

- 3.51 

56.10 

8.18 

3.54 

(V 


7.64 

8.36 

3.53 

1.96 

b) 

5*02 

4.87 

221.7 

57.1 

.1964 

- 2.38 

70.8 

7.91 

(a) 

(bj 


6.83 

7.80 

2.99 

1 

(a) 

b 

3.61 

3.69 


% solution obtained, 
^o solution attoiipted. 




^ta not obtained, »odel destroyed, 
^ot Oeflned- 
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TABLE III.- TEST COMDITIONS POE 45° SWEPT WING 


3 

(a) 

I 

V 

q 

f 

Comments 



214.3 

55.2 

25.3 

Spar locked 

0 

0.1092 

217.2 

56.2 

23.0 

Flutter was "bending- torsion type 

-I.l6 

.1774 



— - 

Low-reduced- frequency flutter not 

-4.098 

-5-88 

.2160 

.2800 


— 

— 

obtained and data for bending- 
torsion flutter were not recorded 


s ie nondimenBloiial in semlchords which are normal to leading 

edge. 
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Figure 2.- Flutter coefficient against Xving location for wood wings I 
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